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We show that the anomaUes in the B — > tttt decay channels, viz., the large branching ratio of 
B — > tt^tt" and large direct CP asymmetry of -B ^ tt+tt", can both be explained, even in the 
framework of existing theoretical models which predict the relative strengths of the tree and the 
strong penguin amplitudes, if there is an abnormally large electroweak penguin contribution. We 
also critically examine the expectation of a similarly large electroweak penguin in the B — > nK 
sector, and show that it can be accomodated, and may even be motivated from an SU(3) flavor 
symmetry between the B — > nn and the B — > nK amplitudes, but is not necessary to explain the 
data. We emphasize that if the experimental numbers on these charmless nonleptonic channels 
survive, they may turn out to be an indirect signal of new physics. Current data is insufficient to 
discriminate betwen flavor-specific and flavor-blind new physics scenarios. 



PACS numbers: 13.25.Hw, 14.40.Nd, 12.15.Hh 
I. INTRODUCTION 

Nonleptonic B meson decays pose a serious theoretical 
challenge to us. The main difficulty is to handle the un- 
certainties of low-energy QCD; to be very precise, to cal- 
culate the long-distance contribution {MiM2\'Heff\B). 
There are at least three models with varied degrees of 
sophistication to tackle the said uncertainty: the conven- 
tional factorization (CF) model jly2i, the QCD-improved 
factorization (QCDF) model 0,1^ and the perturbative 
QCD (PQCD) model These models provide defi- 

nite theoretical predictions for branching ratios (BR) and 
CP asymmetries {Aqp) for nonleptonic charmless B de- 
cay modes. We still do not have comparable numbers 
from other approaches, e. 17., the soft-coUinear effective 
theory (SCET) model Q, though they have predictions 
for some of the B — > tttt amplitudes. 

One of the major obstacles is to calculate the strong 
phase, and hence the CP asymmetry, theoretically. (The 
direct CP asymmetry depends both on the strong phase 
difference and the weak phase difference of the two in- 
terfering amplitudes, and is zero if any one of them van- 
ishes.) The color-transparency argument of Bjorken J3| 
is not expected to hold beyond the color-allowed tree de- 
cays. In the CF model the strong phases are computed 
from the imaginary parts of the respective Wilson coef- 
ficients 0], and no soft final-state interactions are taken 
into account. Among the more sophisticated versions, 
QCDF predicts a rather small strong phase difference 
between the dominant amplitudes. PQCD predictions 
are different as annihilation and exchange topologies are 
given more weightage than in QCDF (and the penguin 
amplitude is also enhanced), and they can generate a 
sizable strong phase. (A comparative study of these two 
models can be found in 4]). This is why the predictions 
for CP asymmetries in QCDF and PQCD are sometimes 
even opposite in sign. Needless to say, the relative impor- 
tance of different amplitudes varies from model to model. 
In particular, strong penguins are much larger in the 



PQCD model due to a dynamical enhancement. There 
are, however, some common trends, which can be justi- 
fied even intuitively. For example, the tree-amplitude in 
B TT+TT" dominates over the strong penguin, whereas 
the situation is opposite for B 7r~K~^. This is en- 
tirely due to the relative importance of the CKM matrix 
elements. 

All these models also agree on the fact that the elec- 
troweak penguins (EWP) are much smaller than the 
strong penguins. That can be seen very easily, by com- 
paring the relevant Wilson coefficients. Intuitively, it 
follows from the smallness of a compared to the strong 
coupling constant ag. In fact, one can safely neglect all 
the electroweak penguins except that mediated by a top 
quark in the loop. However, one must remember that the 
strong phases coming from gluonic penguins (which is 
necessarily / = 0) need not be the same as those coming 
from 7, Z penguins (they do not have a definite isospin). 
Again, this complication does not matter since the EWPs 
are supposed to be very small. There is no known mecha- 
nism which can enhance the EWP amplitude to the level 
of, say, the corresponding gluonic penguins. 

The predictions from the theoretical models are sup- 
posed to be robust when both the daughter mesons are 
light. For example, BRs and Acps have been computed 
for all the B — s- tttt and B — > ttK modes. Such numbers 
mostly agree with each other, if we take into account the 
theoretical uncertainties. This leads us to believe that 
the leading predictions of these models may be taken se- 
riously. These include the relative importance of the tree 
and strong penguin amplitudes, the prediction for strong 
phase differences, et cetera. 

However, as is well known and will be substantiated 
in the next section, some of the experimental numbers 
are at variance with the theoretical predictions, though 
most of them agree. In the B — > tttt sector, there are 
three such discrepancies: (i) Br{B tt'^tt^) is abnor- 
mally large (the largest prediction is from PQCD, which 
is about 5 times less); (ii) Br{B tt+tt^) is slightly 
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on the lower side (the lowest prediction is again from 
PQCD, which is more than 2cr above the data; QCDF 
has a large theoretical uncertainty and taking every- 
thing into account, may just fit the data); and (iii) the 
direct CP asymmetry in B —^ tt+tt" is larger than theo- 
retical prediction, irrespective of the model chosen. One 
must be cautious: the error bars, particularly for the CP 
asymmetry data, are still large, and it may be too early 
to draw any definite conclusions based on them. Still, 
the trend is worth investigating. 

In contrast, the individual BRs and Acps in the B — s- 
ttK sector are more or less in agreement with the theory. 
But when one eliminates the theoretical uncertainties by 
taking suitable ratios of the BRs, it is claimed |9lll(illll| 
that the EWP contribution is abnormally large. It is 
more puzzling since it appears that we understand the 
radiative decay 6 ^ 57 fairly well (though there are ways 
to evade this problem) . We will critically investigate this 
issue also. It will be shown that solutions with vanish- 
ingly small EWP amplitudes can be obtained even with 
a rough SU(3) flavor symmetry between the tree and the 
strong penguin amplitudes of i? — > nn and B — > ttK 
channels. The strong phase difference between these two 
amplitudes for B — > ttK is small (modulo tt), but is dif- 
ferent from the analogous quantity for B — > rnr decays, 
which is exactly what is predicted in PQCD. However, 
a substantially large EWP can be easily accomodated in 
the data. But there is no apparent conflict between the 
B — » ttK channels and the radiative decay b —^ sj. 

There are two ways of analyzing the data. First, one 
can take a particular model of his or her choice, and fit 
the data to get an idea of the relatively poorly known 
parameters. The trouble is that equal justice cannot be 
done to all the experiments. (A well-known example is 
the pre-2003 CKMfitter fit of the B tt+tt" asymmetry 
using the BaBar data only, within the context of QCDF; 
the Belle data was so far away that it could not be fit- 
ted.) The second approach |3, [13 is to analyze the data 
without paying any heed to the models, except, maybe, 
for some basic symmetry ideas like flavor SU(3). As ex- 
pected, the fitted amplitudes and strong phases, all a 
priori unknown, come out to be different from the model 
predictions. 

The question is whether these discrepancies are due to 
the fact that we do not understand the low-energy QCD 
well, or whether there is new physics (NP) beyond the 
Standard Model (SM). If there is NP, that wiU generate 
some more effective four-Fermi operators, leading to new 
contributions in B decays. If these operators are flavor- 
blind, they should mimic the strong penguins, whereas 
the flavor-specific operators should appear more like a 
modiflcation to the tree or the EWP amplitudes. The 
exact structure, of course, will depend on the specific 
model. Unfortunately, due to lesser number of -B — > tttt 
channels, it is not possible to extract an observable which 
will clearly point to some abnormality in the EWP sector, 
and we have to look for an indirect answer. 

The large BR of the tt^tt" channel tells that if un- 



known EWP dynamics is the cause of the enhancement, 
the EWP amplitude should be quite large (we will later 
see that it is almost as large as the color-allowed tree 
amplitude) . It is difficult to have such a large amplitude 
from intuitive arguments. The same may be said for the 
B ttK channels. (That is why we prefer a NP solution 
for the puzzle.) 

This brings us to the justification of our approach. We 
assume that the model predictions arc sensible, as far as 
the tree and strong penguins (and annihilation and ex- 
change diagrams) go, and as a typical example we take 
the PQCD model. Qualitatively same results are found 
for other models too, but there is a reason why we choose 
PQCD. As we will show, a fit to the B ^ tttt data needs 
large penguin contribution, and only PQCD can provide 
such a large enhancement, so that this is a more con- 
servative approach to look for NP effects. Also, a large 
\P/T\ helps the angle 7 to be fitted in the first quadrant, 
in accordance with the standard CKM fit • We treat 
the parameters of the EWP sector as unknowns, and fit 
them from experimental data. This will show us that 
it is the TTTT sector, rather than irK, that contains the 
puzzling EWP behaviour; and one should expect some 
deviation from the SM prediction for _B — > /97, instead of 
B K*^. We will also discuss how a large EWP ampli- 
tude can be correlated between B — > tttt and B — > ttK 
modes. 

However, the main emphasis lies elsewhere. Such an 
analysis is the first step to look for effects coming from 
NP, since a number of NP models mimic EWP dynam- 
ics, by having unequal strength of 6 — s- uud and h ddd 
amplitudes. There are numerous examples, starting from 
models with an extra Z' |14| and/or exotic quarks {e.g., 
vector singlets) to R-parity conserving and violating su- 
persymmetry . We do not go into any of these specific 
models; this is a model-independent study to show that 
such models may have an interesting future. 

We wish to point out that this study implicitly assumes 
the unitarity of a 3 x 3 CKM matrix (thus, allowed ranges 
for the amplitudes for a model with a fourth chiral gen- 
eration should be different). A weaker assumption is no 
NP effect in — B^ mixing. Such an effect modifies 
the phase of the -B" — B^ box, and f3 as measured from 
B — > J/ipKs need not be the true f3. This also weak- 
ens the limits on Vtd- Inclusion of such NP effects in 
— 50 mixing does not invalidate our conclusions; it 
only changes the allowed range of the angle 7 of the Uni- 
tarity Triangle (UT). 

The paper is arranged as follows. In Section II, we 
tabulate the experimental and theoretical inputs for the 
analysis. Section III and IV deal with the B ^ tttt and 
B — !■ t:K data respectively. In Section V, we summarize 
and conclude. 
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II. THEORETICAL AND EXPERIMENTAL 
INPUTS 

Let us first set our notations and conventions. We 
work in the {a, f3, 7) convention of the UT. We use B to 
indicate a flavor-untagged 5° or B", or even a charged 
B meson when there is no chance for confusion. 

The CP asymmetry for f is defined as 

which agrees with the convention of 0, 0, 0| but is op- 
posite to that used by |^ 0| . The same convention (ac- 
cording to the quark content, b or b) is used for neu- 
tral B decays to flavor-specific final states, like tt^K^. 
For a flavor-nonspecific final state / {e.g., tt+tt^), wc 
define I = exp{-i2P){f\n\B0)/(f\n\B°), and a^p = 
(1 - + a'Sp = 2lml/{l + For the tt+tt" 

system, they are related to the Stttt and C(A)7r7r of BaBar 
and Belle by 



Mode 


CP 


QCDF 


PQCD 




9.0 - 12.0 


3.5 - 15.0 


5.5 - 9.0 




19.7 - 21.6 


-4 - 21 


-30 - -16 




0.35 - 0.63 


- 0.07 


0.2 - 0.4 




-42.2 - 45.9 


-100 - 




tt+tt" 


3.0 - 6.8 


4.9 - 6.1 


2.6 - 5.0 




0.0 - 0.1 


0.0 


0.0 




14.0 - 18.0 


3.0 - 24.5 


13.0 - 18.6 




-10.6 - -3.7 


-13 - -6 


12.9 - 21.9 




5.0 - 7.4 


1.5 - 10.0 


8.3 - 8.9 




24.4 - 36.5 


- 7 






9.4 - 12.0 


2.0 - 16.0 


7.6 - 11.0 




-9.2 - -2.6 


-10.6 - 3.0 


10.0 - 17.3 




14.0 - 22.0 


5.5 - 26.0 


13.7 - 19.7 




-1.5 - -1.3 


-1.75 - 


0.6 - 1.5 



TABLE L Predictions for BRs (first row, multiplied by 10^) 
and Acps (multiplied by 100, and in our convention) of 
different nonleptonic modes. The numbers are taken from 
Q,H,IEIS' The uncertainties are not treated equally; in par- 
ticular, QCDF pedictions take into account more sources of 
uncertainty than the other two. The UT angle 7 is fixed to 
be in the first quadrant. For the modes indicated by J, the 
CP asymmetry is measured with a Ks in the final state. 



We also use the following valence quark convention: 

B° = bd, B+ = bu, K° = sd, 

K+ = su,TT+ = ud, TT° = {uu - dd) / V2, (3) 

and for the antiparticles, the quark contents are just re- 
versed, without any additional — sign. This is different 
from those used by the Gronau-Rosner group 0, but 
is consistent with, say, 0. In the convention of 
the following mesons have an extra minus sign in their 
wavefunctions: B^ , K~ , tt^, tt^. This obviously does not 
change the BRs or Acps. However, they change the am- 
plitude sum rules. This will be discussed later in this 
section. 

The experimental data (updated for Winter 2004), 
taken from the HFAG website jlGjj and quoted at 68% 
confidence level (CL), is as follows. 

The BRs, muhiplied by 10'^, are 

Br{Tr+TT-) = 4.6 ±0.4; 

Br{n°TT°) = 1.9 ±0.5; 

Br{TT^TT°) = 5.2 ±0.8; 

Br{K°TT+) = 21.8 ± 1.4; 

Br(i^+^°) = 12.8t}:J; 

BriK+TT-) = 18.2 ±0.8; 

Br{K\°) = n.9tli. (4) 



The CP asymmetries are 

at^C^^ = -0.46 ±0.13; 

a^^ = -S^^ = 0.74 ±0.16; 

Acpin+TT°) = -0.07±0.14; 

AcpiK°Tr+) = -0.02 ±0.06; 

AcpiK+n") = 0.00 ±0.07; 

AcpiK+TT-) = 0.095 ± 0.028; 

Acp{K"7r°) ^ -0.03 ±0.37; 

= ^Ks.- - 0.40l°;27 ± 0.09; 

^Ks-^-^-Sks.- = 0.48t°i8±0.28. (5) 

The theoretical predictions based on different models 
are given in Table 1. The 'large' BRs (i.e., except that of 
B 7r°7r°) more or less agree with each other; there is a 
nontrivial overlap region. Also to be noted that only the 
BR of -B — > TT^TT*^ is significantly off the mark. The pre- 
dictions for Acps differ, mainly because an enhancement 
of penguins and annihilation topologies in the PQCD 
model. However, even taking into account all the the- 
oretical uncertainties, some of the experimental resuls do 
not tally with the predictions. The error bars are large, 
and ultimately the difference might go away, but this is 
high time to be prepared for any indirect signal of NP. 

Let us first discuss the theory of J5 tttt. We fol- 
low the standard convention of writing the amplitudes 
(i.e., putting the top- mediated penguins together with 
up- and charm-mediated penguins by using the unitarity 
relationship VtbV*^^^^ = -VubV^^/^-VcbV*^/^ except that 
the EWPs, both color-allowed and color-suppressed, are 
written separately. This will help us to understand the 
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EWP dynamics of the tttt sector. The amplitudes are 

+i\PEw\e''-'^^-'^^G 



PEw\e''^^^-^^^G, 



(6) 



where = 1 + i?6exp(— 17), and Ri, = |Ati/Ac| = 
0.37 ± 0.04, with A, = V^bV*^. Conventionally, and C 
contain, apart from the tree-level amplitudes, those parts 
of the strong penguin that are proportional to e"*'''. For 
simplicity (and also guided by intuitive arguments) we 
take only the top-mediated EWP amplitudes; they are 
proportional to At, and can be written in terms of A„ 
and Ac using the unitarity argument. Note that we have 
extracted the strong phase associated with the leading 
term, since only the difference of strong phases is a mea- 
surable quantity. The annihilation contributions have 
been neglected, but can be suitably dumped with Tc or 
Pc- To calculate the respective BRs, one has to remember 
that B+ lives longer than B° by a factor of 1.086 ±0.017 

13 . . 

Due to our slightly different valence quark convention, 

the amplitude sum rule looks like 



A{7r+n-) ~ \/2A(7r*'7r°) = %/2A(7r+7r*'). 



(7) 



This convention, of course, yields identical results to that 
of, say, [13 • 

Thus, there are ten free parameters; five amplitudes, 
four independent strong phase differences and the weak 
phase 7. The observables are far less in number: only 
six, namely, three BRs and three CP asymmetries. One 
therefore cannot say in a model-independent way that 
the EWPs are going to be large. However, if one neglects 
the EWPs, the number of free parameters gets reduced 
by four (two amplitudes and two strong phases), so one 
can extract a unique solution l^. We, on the other hand, 
will take the model predictions seriously. For example, 
we use the PQCD prediction = \PclTc\ = 0.23^^ !^^ 
and -41° <8 = 8p-8t< -32°. We also use the CKM 
fit value for 7: 50° < 7 < 72° at 68% CL. This gives us 
three more theoretical inputs. 



Next we come to the B — > -kK modes. One can in- 
voke the flavor SU(3) and relate the B ^ tttt amplitudes 
to the B ttK ones. The symmetry is supposed to 
be broken only by the corresponding CKM factors and 
the ratio of the meson decay constants ]kI j-n- This in 
turn implies that one takes the chiral enhancement fac- 
tors coming from {V — A)®{y + A) type operators (O5- 
Os in the standard literature) to be the same in both 
cases. Let us follow a different approach. We start with 
the assumption that the amplitudes are independent of 
each other; later we will see whether there is any trace 
of the flavor symmetry. If one can relate the amplitudes 
by SU(3), that will act only as a further boost to our 
analysis. In fact, the PQCD relationship between the 
amplitudes and Pc, and the analogous one for 6 — s- s 
transitions, \T'JP'^\ = 0.201 ± 0.037, is consistent with 
the SU(3) expectation 



T' 

^ C 

P' 



Tc-^— = 0.28Tc; 

J-K Vud 

p Ik_Vc^ = -5.3Pc, 



(8) 



where primes indicate b ^ s transitions. Note that 
PQCD also predicts the strong phase difference between 
Pc and T'^ to be about 156° |6j]. We, however, wiU not use 
that as an input; rather, we will show that the allowed 
solutions, obtained just from fitting the data, are entirely 
consistent with the PQCD prediction. 

We use the Wolfenstein representation of the CKM ma- 
trix, and define A = - P/2) w 0.230, and 



Rb{d 



Vus Vcd 
Vud Vcs 



-X^Rb. 



(9) 



Remember that Vcd is negative, hence the minus sign. 
We assume that there is no effect of NP or unknown SM 
dynamics in B^ — B'^ mixing, and the phase from the B^ — 
P° box is sin(2/3), which has been directly measured as 
well as fitted. We use the CKMfitter fit result sin(2/3) = 
0.739 ± 0.048. It was discussed earlier what happens if 
one relaxes this assumption and allows NP effects in the 
mixing. 

The amplitudes can be written in an analogous way to 
that of the B — > tttt ones. There is, howver, one impor- 
tant difference. The modes B+ K^7r+ and B K^n" 
have no tree-level contributions. The strong penguin, 
which after rearrangement was absorbed into Tc, now 
gives a different amplitude P^t = P^ — Pi with a differ- 
ent strong phase S'^jj, (we use primes to indicate 6 — > s 
transitions). 
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A{B K-Tr+) 
A{B- ROtt-] 



.\P'\ — IP' \p-i7pi(S'cT-S'p) J- ^\pC '\„i{S'Ec-S'p)Q' 

A ^ 

|P^'| + |P4|e-'^e'(''cT-'5i.) _ i|pC^'|ei('5Ec-'5p)^' + 
A ^ 



V2A{B- K--K°) = |r; 



V2A{B^K0Tr°) = 



(10) 



In the above expressions, G' = l + i?^e~'^, and C is the 
same as C (in B mr) except that there is no strong 
penguin (this difference is non-neghgible) . In the last 
two amphtudes, there is a part where the tt is emitted 
from the B meson; this, expectedly, is very similar in 
structure to that of P ^ tttt. In these terms, the relative 

sign between A and A terms is positive whereas in the 
analogous expression for B nn, that is negative; this 
is because Vcd is negative. Pew is, of course, the same 
as that in i? ^ tttt. Thus, there are five new amplitudes 
and four strong phases. 

III. ANALYSIS: B ^ mv 

The main emphasis of the B tttt analysis is to ex- 
tract, unambiguously, a signature for large EWP ampli- 
tudes. In the SM, we expect Tc to be the largest ampli- 
tude, while Pc and C are smaller than Tc- In fact, models 
like PQCD predict \Pc/Tc\ ^ 0.25, while the predictions 
from other models are generally below 10%. There is 
no such prediction available for C, but one can expect 
\C/Tc\ < 1. The EWP amplitudes are expected to be 
highly suppressed compared to these. The color suppres- 
sion further diminishes Pew compared to Pew- 

For those readers who do not want to go through the 
nitty-grittics of the analysis, let us state our conclusions 
right here. We will see that the amplitudes Tc and C turn 
out to be as expected {Pc is related to Tc), and though 
the color-allowed EWP Pew may happen to be large, 
there arc a number of solutions (we have more unknowns 
than equations, so there is a large number of possible 
solutions) where this turns out to be small, almost as 
small as expected from the theory. Later, we will see 
that only those small values are allowed by the B 
nK data. The surprise comes from the color-suppressed 
EWP amplitude Pew''' turns out to be large, larger 
than Pc, and almost as large as Tc! There arc other 
interesting features of the analysis; we discuss them at 
the end of this section. 

We start by treating all the amplitudes as free param- 
eters, with the exception of Pc- The latter is varied over 



its entire range allowed by the theoretical uncertainties 
of \Pc/Tc\. The other four amplitudes are varied over the 
range 0-6 x 10~^. We will see later that this is a rather 
conservative range. All the strong phase differences are 
varied over the range 0-2^, except Sp — dx, which is again 
constrained by the theory. The weak angle 7 is varied be- 
tween 50° and 72°. 

From the solutions that satisfy the experimental con- 
straints, we get the following ranges for the amplitudes: 

2.5 < Tc X 10« < 4.1; 

0.4 < Pc X 10^ < 1.2; 

< C X 10« < 3.8; 

< Pew X 10^ < 2.5; 

1.0 < P^w X 10* < 3.7. (11) 

Note that apart from the last line, none of these are 
unexpected. One may object that C is supposed to be 
smaller than Tc, if we take that constraint into account 
nothing changes except the upper limit of C; it becomes 
3.5 X 10~*. However, there is a concenration of solution 
points in the C — Pew plane near the origin (see Fig. 
1). This vindicates the conjecture that both these ampli- 
tudes may be small. Later we will see from the B ttK 
channels that the upper limit of Pew is smaller than the 
range predicted here by about 40%. 

Why Pew i^ large? The main culprit is the B 
7r°7r° BR. This, along with the facts that the other two 
BRs arc more or less in the expected ballpark, while the 
direct CP asymmetry in i? ^ tt+tt" is large (i.e., there 
must be a significant interfering amplitude) generates a 
large Pew- ^iS- 2 to have an idea of the solution 
points for P^W' i^ ^^^^t, if we claim that Pew is small, 
the range for P^^y becomes 2-3x10"^ This is a large 
amplitude by any standard, and cannot be accomodated 
in any theoretical model. 

The strong phase differences do not show any inter- 
esting pattern, except, again, for 6ec — ^t- Varied be- 
tween and 27r, solutions come out only for the range 
174° < Sec — St < 190°. This shows a large destructive 
interference between Pew ^^'^ -^c) which is necessary to 
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FIG. 1: C versus Pew for B mv. Note the concentration 
of solutions near the origin, and the excluded region for Pew 
for low C. The amplitudes have been multiplied by 10*. 




FIG. 3: Direct and mixing-induced CP asymmeries in _B ^ 
7r"'"7r~. The circular arc that bounds the parameter space at 
the upper left hand corner corresponds to la^JirP + |i™7rP = 1- 



reproduce the BRs of i? ^ tt+tt and i?+ tt+tt" in 
the right ballpark. 




FIG. 2: Pew versus Pew foi" ^ Note the large nonzero 

allowed values for the latter, and the concentration of solu- 
tions for small Pew- The amplitudes have been multiplied 
by 10«. 

Fig. 3 shows the allowed parameter space of a^^ ver- 
sus o™^ ; this shows a concentration of solutions near the 
physical boundary |a^^p -I- |a™^p < 1. 

Let us emphasize again at this point that this treat- 
ment is valid if and only if the new physics, if any, con- 
tributes to the decay but not to the B^—B^ mixing. Most 
of the new physics models are of this type; however, there 
are exceptions. For example, four-generation fermion 
models contribute mostly in mixing. Supersymmetry 
with R-parity violation may contribute to only mixing, 
only decay, or both, depending upon the nonzero cou- 
plings chosen. Generically, if NP contributes in B^ — B^ 
mixing, the following modifications have to be taken into 
account (an example with R-parity violation has been 
exhaustively discussed in 18] ): 

• A Aid is affected by NP. Consequently, Vtd, which 



is determined from the measurement of AM^, be- 
comes a free parameter, except for constraints com- 
ing from the unitarity of the CKM matrix. Obvi- 
ously, for models with more than three quark gener- 
ations, this is a very loose constraint. Fortunately, 
the precise value of Vtd is not that important in our 
analysis. 

• What is more important is the value of sin(2/3). 
There is a SM value of /?, which let us call Psm- 
We have explicitly assumed that the CP asymme- 
try in, say, B J/^Ks, measures /3sm- If there 
is NP in mixing, it will measure f3 which should be 
different from (3s M- Note that this happens even 
if the NP mixing amplitude is real. The value of /3 
is important in the analysis of CP asymmetries in 
the B TT+TT" channel. To get a feehng of Psm, 
one should perform the usual CKM fit for the ver- 
tex (p, 77) of the UT without the AMd constraint. 
Again, this is not a foolproof prescription; one im- 
plicitly assumes that NP does not contribute to, 
say, e'/e. 

Even considering all this, it is hard to see how a NP 
that contributes only to mixing but not to decay can 
explain the B ^ tttt results. It can affect the fit for 7, 
but cannot generate a large BR for i? — tt^tt*^ or a large 
direct CP asymmetry for B —^ Tr~^n~ . 



IV. ANALYSIS: B ~> ttK 

The analysis for B — > ttK follows that of S ^ tttt. 
There are five new amplitudes, apart from Pew, and 
hence four new strong phase differences. We scan this 
ten-dimensional parameter space for possible solutions. 
Note that i?^ is a small number and hence for all practical 
purpose G' can be approximated by unity (but for our 
analysis we keep it as it is). Also note that these ten 



parameters arc not really all free; the allowed range for 
Pew, a-s well as the weak phase 7, are determined from 
the B ^ nn data, and we use the PQCD result IT^'/P^I = 
0.201 ±0.037. We also equate the strong phases of and 
P' 

We obtain the following ranges for the amplitudes: 



0.2 < 


X 10^ 


< 1.15 


0.8 < 


P^ X 10^ 


<4.8; 


< 


P4 X 108 


< 3.2; 


< 


C X 10^ 


< 3.8; 


< 


Pew X 10« 


< 1.55 


< 


PEw' X 10« 


< 6.3. 



(12) 



There are two points to note right here. First, C could 
in principle have a bigger range. We constrained it to 
be smaller than C of S ^ tttt modes. Second, both the 
electroweak penguin amplitudes can be close to zero; we 
refer to Figure 4 for a feeling of the solutions. Thus, 
there is no apparent EWP anomaly in the data. We have 
explicitly checked that solutions are possible even if one 
puts by hand both the EWP amplitudes equal to zero. 

All the strong phase differences have been varied be- 
tween and 27r. Solutions are allowed for the entire 
range, unlike B tttt, where wc found a rather con- 
strained fit 174° < 5ec - St < 190°. However, from 
Figures 5 and 6, there appears to be an interesting cor- 
relation between the amplitudes and the strong phases. 




FIG. 4: Allowed solutions for color-allowed and color- 
suppressed EWP amplitudes, showing that there is no appar- 
ent EWP anomaly in B — » nK data. The amplitudes have 
been multiplied by 10*. 



From Figure 5, it appears that small Pew' solutions 
are admitted only if is large: |T^| > 0.55 x 10"^. For 
smaller values of T^', Pew ™ust be large. This is intu- 
itively clear: after all, one of them should bo largo enough 
to generate the required BRs. One could have done this 
taking and P^. to be free parameters; but we wish to 
follow the theoretical prediction, which is also consistent 
with an SU(3) flavor symmetry, as far as possible. 




0.2 0.3 0.4 0.5 0.6.^, 0.7 0.8 _.. - 1.1 1.2 



FIG. 5: versus P^w' for B ttK. Note how small P^w' 

solutions necessarily corospond to large which is consistent 
with SU(3) relationship. The amplitudes have been multiplied 
by 10^ 
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FIG. 6: Allowed solutions for the strong phase difference 
5p — S't as a function of T^, multiplied by 10®. 



Figure 6 shows the solutions for the strong phase dif- 



feronco S'p — 8'r, 



, function of T'. Note the band aroimd 



the theoretical prediction of 156° from the PQCD model, 
which appears only for large and hence small EWPs. 
This shows that there is nothing to worry about the EWP 
sector; but one can accomodate large values, which may 
come from NP models. 

Thus, the SU(3) flavor symmetry is there, as far as 



.(') pC) 



and Pew concerned. 



the magnitudes of Tc , Pc 
However, the strong phase differences are not necessarily 
small. If one demands a correlation between both the 
amplitudes and the strong phases for B tttt and B 
ttK, one must invoke large EWP contribution. But such 
a contribution is there for B — > tttt; there is no a priori 
reason why it should not appear for B — > ttK, and thus 
help to salvage the flavor symmetry. 
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V. SUMMARY AND CONCLUSIONS 

We have juggled with a number of amphtudes and 
strong phases. The equations being horrendously cou- 
pled, we do not even envisage to attempt an analytic 
solution. Those solutions necessarily involve some sim- 
plification and run the risk of missing the point. Rather, 
we solved these equations numerically, and the equations 
being more in number than constraints (and remember 
that these constraints often have large error bars), found 
a range of allowed values for the free parameters. It is 
indeed interesting that one can form some idea about the 
nature of the solutions. 

We assume that the theoretical models can catch 
the essence of these nonleptonic transitions successfully. 
There are points where they differ; for example, the 
way how one takes care of end-point suppresions, higher- 
twist corrections, or annihilation topologies. Fortunately, 
such subtleties are not really crucial for our analysis. 
What is important is the relative strength of the strong 
penguin amplitudes with respect to the tree amplitude. 
PQCD, by virtue of a dynamical enhancement mecha- 
nism, predicts a larger \Pc/Tc\ than other models like CF 
or QCDF. Since our strategy is to extract the minimal 
set that docs not fit in these model predictions, we take 
PQCD as our working model. If one takes, for example, 
QCDF, one needs a larger color-suppressed EWP to ex- 
plain the S — > TTTT data. PQCD also helps to fit 7 in the 
first quadrant, which is compatible with the CKM fit. 

We found that the explanation of the B tttt data 
needs a large color-suppressed EWP contribution (the 
color-allowed EWP may be large, but need not be). This 
is entirely due to the large BR oi B ^ ■n^n^ . If this 
data would have been in the expected ballpark, the large 
direct CP asymmetry in i? — > 7r"'"7r~ could be explained 
by a much smaller EWP contribution. Such a large EWP 
contribution also means that there should be something 
worth watching in i? — > p7. 

The situation is not so drastic for B ttK. Even if 
one takes the SU(3) predictions of and seriously, 
there are solutions with vanishingly small EWP ampli- 
tudes. Of course, there is a price to pay: the other two 
amplitudes C and become essentially free, except for 
the constraints that \C\ < \C\ and S'p — S'ljrp. We have 
tried to see what happens to the fit if all amplitudes 
and strong phases are taken from some model (we used 
CF, because that predicts all the amplitudes and also 
the strong phases coming from the imaginary parts of 
the respective Wilson coefficients), and there one indeed 
requires a sizable EWP contribution! 

So, what is the lesson from the B ttK analysis? 



Probably that it is yet early to say that there is any defi- 
nite evidence for large EWP amplitudes, but they can be 
easily accomodated, as is evident from our figures. All 
the penguins are expected to be enhanced by about a fac- 
tor of 5 when going from B — > tttt to B —^ ■nK . This trend 
can be accomodated easily for the color-allowed EWPs, 
but appears to be very hard for the color-suppressed ones. 
Thus, even if the color-allowed EWPs can be explained 
within the present models, the color-suppressed one poses 
a serious challenge. 

Is the situation so bizarre because we do not under- 
stand the dynamics of penguins, or to be more general, 
that of low-energy QCD? No one knows, but from intu- 
itive arguments it is very hard to establish such a large 
EWP amplitude compared to the strong penguins. One 
should not also forget the possibility that the experimen- 
tal data may change over the years, though such an as- 
sumption does not help the analysis. 

Is the solution to be found in physics beyond the SM? 
Again, it is early days, but one may be hopeful, more so 
because this may provide an indirect evidence for such 
NP before the Large Hadron Collider (LHC) comes into 
operation. NP models generate new operators in the ef- 
fective Hamiltonian, or modify the Wilson coefficients for 
the existing operators (an example is the Universal Extra 
Dimension model of Appelquist et al ^19|). If the inter- 
actions are flavor-specific, they will appear as a modifi- 
cation to the tree or the EWP amplitudes. These mod- 
els include exotic quarks, new Z' gauge bosons, R-parity 
conserving and violating supersymmetry and so on. 

A couple of cautions here. First, one needs the exact 
structure of a model to be more specific than the present 
analysis. What one can do is to eliminate models: for 
example, supersymmetry with flavor alignment will not 
generate such a large FCNC amplitude. The exact weigh- 
tage of new operators will also depend on the particular 
model chosen. Second, this analysis has to be slightly 
modified if there is new physics in B'^ — B^ mixing. This 
has been discussed in Section 3. 

Where should one get confirmatory signals for such NP 
effects? Again, that depends on the specific model. We 
have already mentioned B pj; another place would be 
the semileptonic b ^ d decays. For nonleptonic decays, 
one should carefully analyse the B — > /jtt, B — > pp, Bg — > 
ttK, and Bg KK channels. An analysis of angular 
correlations in B ^ pp should be useful. 
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